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A correlation between the highest energy Cosmic Rays (above ~ 60 EeV) and the 

■ distribution of active galactic nuclei (AGN) gives rise to a prediction of neutrino 
^ , production in the same sources. In this paper, we present a detailed AGN model, 

\ predicting neutrino production near the foot of the jet, where the photon fields from 

' the disk and synchrotron radiation from the jet itself create high optical depths for 

proton-photon interactions. The protons escape from later shocks where the emis- 
\ sion region is optically thin for proton-photon interactions. Consequently, Cosmic 

■ Rays are predicted to come from FR-I galaxies, independent of the orientation of 
. the source. Neutrinos, on the other hand, are only observable from sources directing 

their jet towards Earth, i.e. flat spectrum radio sources and in particular BL Lac 
type objects, due to the strongly boosted neutrino emission. 



Abstract 



Key words: UHECRs, radio galaxies, neutrinos, jet structure 
PACS: 98.70.Sa, 98.54.Gr, 14.60.Lm 



Corresponding author. Contact: julia.becker@physics.gu.se, phone: +46-31-7723190 



Preprint submitted to Astroparticle Physics 18 December 2008 



1 The underlying AGN model 



The evidence for a correlation between the arrival directions of ultra high 
ene rgy Cosmic Rays (UHECRs) wit h active galactic nuclei, as reported by 
the lAuger Collaboration! (|2007|. 120081) . supports many long standing expecta- 
tions (jCinzburg fc Syrovatskiil . ll964l ). The active galactic nuclei with the most 
detailed availa ble theory to actually acce l erate protons to beyond 10^° eV are 
radio galaxies flBiermann Sz Strittmatterl . 119871 ) . In all radio galaxies the feed- 
ing of outer radio emitting regions is done via a relativistic jet emanating 
from near a black hole. Shock waves in such jets can accelerate particles just 
as shocks in the Solar wind Shocks in the jet emanating from near the 

black hole start around a few thousand gravitational radii, as it was shown with 
detailed spectral fits of the entire electromagnetic spectrum, including the spa- 
tial structure at the wavelengt hs where it is known . The shocks end as strong 
shocks at kpc or further out ( Markoff et al. . 200 ll . 2005 ). It was also shown 
that when particles get accelerat ed at the first shock, proton-photon interac- 
tions limit their maximal energy f Kriills . 1992 : Nellen et al.l . 1993 : Mannheim . 
1995h . 



The correlation with the distribution of active galactic nuclei cla imed by the 



Auger collaboration has been disputed by the HiRes collaboration flAbbasi et al. 



2008al ). However, it is not clear at this point, whether both data sets use the 
same lower energy cutoff with the same sh arpness. This is important, as the 
MHD simulations of cosmic magnetic fields (IRyu et al.l . l2008l : iDas et al.l . |2008| ) 
show that scattering of ultra high energy particles rapidly increases with lower 
energy even near 60 EeV. Thus, with even a slight mismatch between the two 
data sets the statistics could be very skewed. Using a complete sample of radio 
galaxy sources and their predicte d properties as UHECR sources, these statis- 



tics will be e xplor ed elsewhere ( ICaramete et al.l . l2008l : ICurutiu et al.l . 12008 



Dutan et all . l2008f ) 



1.1 FR- 1 galaxies and UHECRs 



R adio galaxies with exten ded radio jets were classified into two categories 
by lFanaroff fc Rileyl (119741 ): A population of high luminosity shows radio lobes 
at the outer edge of the jet, at kpc scales from the core, FanarojJ Riley II 
objects, short FR-II. The lower luminosity population, on the other hand, 
has radio knots distributed along the jet, the first knot being as close as 
~ 3000 Schwarzschild radii from the central core, Fanaroff Riley I objects, 



^ For a first discu ssion of the Solar wind, see ( Biermann . 195ll ). for shock acceler- 
ation in AGN see (jBednarz Qstrowskil . Il998l . e.g.). 
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short FR-I. Both radio lobes and knots show non-thermal radio spectra, aris- 
ing from electron acceleration at a shock front as first theoretically described 



by iFermil (119491 . Il954l ). In analogy to processes in the creation of Galactic 



Cosmic Rays, protons are believed to be accelerated at those shock fronts 
in the same way as electrons. In particular, oblique shocks can be very effi- 
cient in particle acceleration, and use electric fields in shock-drift acceleration 
due to the Lorentz tran s formation of t he magnetic f i elds in the proper frame , 
see (Hoffmann fc Tellerl. Il950l: iParkerl. Il958l: iJokipiil. 119871: IVolk &: Biermannl . 



19881 : iBiermannl . Il993l : iMeh fc Biermannl . l2006l : iMeh et all . l2008f ) 



1.1.1 Unified model of FR-I galaxies and BL Lac objects 



In this paper, we consider FR-I galaxies as the so urces of the UHECRs poten- 

tially observed by Auger. It is discussed by others (Tascau ^ 2003[ ^ Cuoco fc Hannestad 



2008 



2008 



Halzen k O'Murchadhal . l2008l : iKoers fc Tinvakovl . l2008l : iKachefriess et al. 



that the nearby FR-I galaxy NGC5128, Centaurus A (Gen A in the fol- 
lowing), is a good candidate to be responsible for a large fractio n of the corre- 

lated events above 60 EeV. M 87 is another closeby candidate, see lBiermann fc Strittmatter 
( 119871 ). which cannot contribute to a possible Auger correlation, since it is 
barely in Auger's field of view. The large number of more distant FR-I galax- 
ies provide good candidates for the total Gosmic Ray flux above the ankle. 
Here, we discuss the morphology of FR-I type galaxies and how these can 
accelerate particles to the highest energies. 



Figure 1 presents a schematic view of the model of FR-I galaxies that we use 
in this paper. On the x-axis, the rotational-symmetric part of the AGN is 
shown, while the y-axis represents the axis of rotational symmetry along the 
AGN jet. Both axes have logarithmic units. When the AGN jet is pointed 
towards Earth, the FR - I typ e galaxy is viewed as a BL Lac type object 
( lUrry fc Padovanlll994l . Il995l ). showing fiat radio spectra, with an unresolved 
jet structure. When the AGN jet is viewed at an angle, the jet structure with 
radio knots distributed along the jet can be observed. In contrast to the more 
radio-luminous FR-II galaxies, FR-I type objects typically lack the observa- 
tion of optical disks. While some of the FR-I type galaxies, such as M 87, 
clearly lack luminous accretion diskj^ and tori, many objects in this class 
may have accretion disks hidden behind the tor us, provided that the torus is 



closed around the jet as indicated in Fig. 1, see (IFalcke et al.l . Il995al ). In this 



case, the jet-disk symbiosis model holds also for FR-I galaxies and the disk 
power scales with the radio luminosity. 



^ They are likely to have radiatively inefficient accretion disks, typically for low 
accretion rates. 
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BL Lac 




Fig. 1. Schematic figure of the class of FR-I galaxies with double- logarithmic scales. 
When the jet is pointed directly towards Earth, FR-I galaxies are classified as BL 
Lac objects. In this model, the torus will hide the accretion disk from view. In the 
case of a FR-I galaxy without torus, no radiatively efficient disk may be present. 

1.1.2 Magnetic fields and shock structure in FR-I galaxies 

The dependence of the magnetic field in these jets along the jet axis Zj is 
near B ~ at large distances. However, in the inner region, the radial 
dependence is not certain. Here, we investigate the radial dependence of the 
magnetic field in order to determine the Cosmic Rays' maximum energy along 
the jet and in particular, at the innermost part of the jet. 



Observations (IBridle fc Perleyl . 11984 e.g.) suggest that the magnetic field runs 



as z~'^ at first, since the radio polarization observation show that the mag- 
netic field is paral lel to the jet. By analogy to the solutions of divB = in 
a magnetic wind (jParkerl . Il958l ). the magnetic field locally shows a parallel 



component of ~ Zj ^. Further out, the magnetic field observations suggest 
that, just as in a wind, the magnetic field becomes perpendicular to the flow 



direction, and so B 



It is obvious that in a smooth wind. 



any compo- 
nent decreasing with z~^ will ultimately win over a comp onent running as z'"^ . 



As the equation o f state is almost certainly relativistic (IFalcke fc Biermanru . 
19951 : iFalcke et al.l . Il995bl . e.g.), the total pressure P depends on the density 
as P ~ P^^"^ 1 while in a conical simple jet the density p ~ giving in near 
equipartition (the magnetic field pressure running with the total pressure) 
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then B 



P(r) 



-4/3 



Howev er, repeated shock waves will reheat the material (see e.g. ISanders 



( 1983 ) , as w e ll as Mach's original work in the 19th century, iMach fc Wentzel 



1884ll885h : lMachl fll898h ). and so we will assume that the Mach-number of 



the flow repeatedly comes back to the same value, while the jet flow velocity 
will stay approximately constant. Therefore going from crest to crest 



P 



and so as a consequence 



(2) 



This is consistent with the concept that the jet stays approximately conical. 
This argument is independent of the orientation of the magnetic field, and so 
the radio polarization observations are not in contradiction, but need then an 
interpretation as arising from highly oblique shocks, which emphasize magnetic 
field components parallel to the shock surface. Highly oblique shocks are only 
possible for high Mach-numbers, which again is consistent. 



stream lines 
Stationary Shoci< Front 
Jet boundary 

stioci< angie towards jet boundary 



moving stiocl< region 
~ 30 r ~ 3000 r 



/ 



~ 3000 r 



stationary stiocl< region 



Fig. 2. Proposed conical shock structure in AGN jet. 

As one check let us consider a system of repeated conical shocks, and ig- 
nore for_^implicity the inner Mach disks as shown schematically in Fig. 2 (see 
also I Sanders! (119831 )). Then we can see, that given a specific Mach-number 
at the initial flow formation there will be highly oblique shock waves, which 
will repeatedly reflect on the conical boundaries, and so produce a self-similar 
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pattern as long as the Mach-number keeps returning to near its initial value. 
A self-similar repeated structure with an ever increasing inner scale will re- 
sult. This is precisely what is seen in jet structur e at vastly discrepant s patial 
resolutions, like in the radio galaxy NGC6251 (IBridle &: Perleyl (119841 ). and 
earlier papers). Further examples for suc h a sh ock structure are recent obser- 



vations of BL Lacertae (IMarscher et al.l . |2008| ) and the BL Lac type object 



S5 1803+784 flBritzen et al.l . 120081 ). 



For those oblique, stationary shocks, it is likely that a mixture of sub- and 
superluminal shocks is presenif^. The proton spectra look very different when 
comparing the two cases. While superluminal spectra have maximum energies 
of around 10^ GeV, subluminal shocks ca n reach energies up to the highest 
energies, i.e. -Emax ~ 10^^ eV, as shown by iMeli et al.l (120081 ). 



Another consequence is that the magnetic field can be sufficiently high even 
far out to confine particles at the energies of UHECRs. As the magnetic field 
is anchored to the inner accretion rate, or its residual electric currents from an 
earlier accretion event, any magnetic field that decreases much faster than 
on the way out, will have difficulty to confine pa r ticles during their accelera- 
tion t o the highest energies obse rved (IBlandfordl . 1 19761 : iBlandford fc Znajekl . 
19771 : IBlandford &: Konigll . Il979l ). Only a magnetic field which runs overall as 
B 

structure 



Zj allows the magnetic field to be relatively high far out along the jet 



We conclude that the magnetic field then runs as B 



approximately. 



1.1.3 Cosmic Ray acceleration in FR- 1 galaxies 



Normally when considering whe re UHECRs can be accelerated, the spatial 
limit, or Hillas-limit, is invoked (iHillasl . 1 19841 ). Using radio observations this 
suggests that radio hot spots of powerful Fanaroff- Riley II radio galaxies are 
very good bets: they are usually modelled as shocks (iMeisenheimer et al 



1989), which can be shown to accelerate particles to near 10^^ eV (IBiermann fc Strittmatter 
1983). 



When using the information on the magnetic fields inferred from radio jets, 
their r adio knots, and their hot spo t s, one typically finds lO"'' Gauss on kpc 
scales ( iMileyl . Il980l : iBridle fc Perleyl . 1 19841 ). On the other hand, the maximum 
magnetic field close to the black hole, so on scales such as a few Schwar z schild 
radii, is g i ven by 10^ Gauss (Mr^/IO'^ M^)"^/^ (IShakura fc Sunvaevl . Il973l : 



Blandfordl . Il976l : iMassi fc Kaufman Bernadol . |2008| ) . Therefore on radial scales 



^ The value of the angle between the magnetic field and the shock front normal 
determines whether a transformation into the Hoffmann- Teller frame (E = 0) is 
possible (subluminal) or not (superluminal). 
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over a factor of 10^ the magnetic field decreases by just this fa ctor, quite indica- 
tive of a behavior. This was also used quite successfully by lBlandford fc Konigl 
(119791 ) ■ Therefore, we conclude that such a radial dependence of the overall 
magnetic field is well justified. 

When estimating what the magnetic field might be in a jet, we have used the 
curre nt accretion rate to i ndicate an estimate, conn ecting it to the accretion 
disk (IFalcke fc Biermannl . Il995l : iFalcke et al.l . Il995bl ) . This then gives usually 
a rather weak field in all cases when the overall activity is low, such as believed 
in BL Lac objects. And a weak field fails the Hillas test for UHECRs. 



However, an alternative proposed by iBlandford fc Znajekl (119771 ) is that the 
magnetic field is still quite high, from a prior accretion episode, and then the 
jet can be driven by a spin-down of the black hole. In this context the magnetic 
fields are higher, and even lo w power sources such as FR-I ra dio galaxies are 
possible sources of UHECRs (IDutan fc Biermannl . l2005l . 120081 ) . 



In addition, the Lovelace limit (ILovelacd . Il976l ) shows that the Poynting flux 



a lower limit to the energy flux in a jet, is connected to the maximal en- 
ergy of a particle confined in the jet by Lj^t = 10^'^erg/s (E'max/lO^^eV)^. 
Therefore a jet such as Cen A , estimated to carry probably around 10^^ erg/ s 



( IWhysong &: Antonuccil . 120031 ) cannot possibly accelerate particles - especially 
protons - to 10^^ eV, at most it would seem, that 10^^ eV is possible. This limit 
can be exceeded by three arguments This limit can be exceeded through three 
arguments: (a) The particles might be heavier than hydrogen. I n this case, 



(2005 


); 


Ave et al. ( 


2005) 



tions and this argument does n o t app ly here, (b) Secondly and more im- 



portantly, iGallant fc Achterbergl (119991 ) show that the Lorentz factor of the 
shock in the local upstream frame enters squared, (c) Further, the jets might 
be intermittent, as strongly demons t rated by the radio galaxy Hercules A 



(iGizani fc Leahvl . 120031 : iNulsen et al.l . 120051 ). 



Therefore there is no a priori difficulty for FR-I radio galaxies to accelerate 
protons to near 10^° eV. With a magnetic field decreasing as zj^ ^ the Hillas 
spatial limit will give the same maximal proton energy at all radii. 



1.2 Optical depth 



We discuss the interaction of ultra high energy cosmic rays with three different 
targets along the AGN jet. We calculate the optical depth of interactions with 
the photon fields of the disk and the knots' synchrotron field. As a third 
possible target, we consider the interaction of UHECRs where the jet meets 
the torus (see Fig. 1). 
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'1) Proton interactions with disk photons 

The optical depth for proton-photon interactions is given by the ratio of 
the length I k, Zj ■ 9 and mean free path of the protons, Ap^^.^j^, in the jet 
in the disk's photon field n.^^^^^: 



I 



r, 



P7disk 



A 



(3) 



■P7disk 



Here, = 900 yubarn is the total cross sectio n for the production of the 



Delta -resonance in proton-photon interactions (IBiermann fc Strittmatter 



1987f ). 



To check on a realistic value for Tp^^^^^ we estimate the interaction prob- 
ability as follows: 

In an active galactic nucleus the accretion disk will produce a radiation 
field near to the nucleus. The accretion disk luminosity at full efficiency 
is given by 



-^^disk ~ esdd ■ 10 erg/s 



(4) 



where eEdd < 1 is the accretion rate relative to the maximum, the Ed- 
dington rate. 

At some distance zj along the jet, starting at about 3000 r^, where = 
1.5 ■ 10^^ MbhI (lO ^MfTi) cm, the gravit a tiona l radius, we find stationary 
shock waves, see (IMarkoff et al.l . l200ll . |2005| ) and also iMarscher et al. 
(120081). Moving shock waves ar e expected to be presen t between ~ 10 — 
3000 as discussed in detail in lMarscher et al.l (120081 ). These Shockwaves 
accelerate particles, and these particles, say protons, interact with the 
radiation field. The photon density is then given by 



n. 



disk 



7disk 



Ait Zj c ■ hv 



(5) 



where hv is the typical photon energy, for which we adopt 20 eV, which 
is /«/ = 3 ■ 10"^^ erg. Using 3000 as a reference radius, this expression 
can be rewritten as 



n 



4 ■ 10^^ eEdd 



7disk 



disk 



^3 



cm 



(6) 



The optical depth across the jet, with the length across the jet as 9 Zj, is 
the given by 



'P7disk 



0.2 ■ eEdd ■ 



disk 



0.1 



3000 r„ 



10^4 erg/s 



(7) 
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decreasing linearly outwards. Since FR-I galaxies are not as efficient ra- 
diators as FR-II galaxies (see Section 1.1), the Eddington rate in the best 
case is Sedd ~ 0.1. Those FR-I sources with the torus close around the 
jet, the optical depth is therefore ~ 2%. For other FR-I type objects, like 
M 87, where the accretion rate has decreased, the d isk can be come very 



faint . This means changing from a radiative disk, (jShakura fc Sunyaev 



19731: iNovikov fc Thornd . 1 19731 . e.g.), to a radiatively inefficient disk. 



(INarayan fc Yil . ll995l . e.g.), and the optical depth will essentially be equal 



to zero. 

(2) Proton interactions with the synchrotron photon field in the jet 

If protons interact with synchrotron photons in the same knot, the optical 
depth is given as t^^^^^^^^ ~ Zj ■ 9 ■ nl^fl^^^ ■ fp-y. The particle density of 
synchrotron photons, nl^^^^^^^^, is given in the rest frame of the plasma 
and needs to be transformed into the observer's frame, n°^^ . Even if 

fsynch 

the shock is standing, the plasma is streaming with relativistic velocities 
along the jet, and so is the synchrotron photon field. The photon density 
observed at Earth is thus modified from th e field that the interacting 



protons 'see' in the plasma's rest frame, see (IRybicki fc Lightmaru . 1 19791 . 
e.g.): 

n:-* =r-i-nf . (8) 

/synch (syncn ^ ' 

Here, F is the boost factor of the streaming plasma relative to Earth. The 
photon density in the plasma's frame can be determined by assuming that 
the luminosity of a knot along the jet is a fraction eknot ~ 0.1 of the total 
synchrotron luminosity, L^ynch ~ 10"^° erg/s: 

^obs ^ ^knot ■ -^synch /gx 

7synch ^ i-Jl . z'j ■ €■ {hv) ' 

The frequency of synchrotron photons is z/ ~ 1 GHz. Hence, the optical 
depth is 



Tsynch 0.9 



•io\ / e 



f Cknot\ 



F / vo.i ; V 0.1 / 

-^synch \ / Zj 



" ' (10) 



\^ 1040 erg/s/ \?>mQrg) VlGHz 

Thus, for a relativistically streaming plasma of F ~ 10, optical depths 
around unity are expected. At the foot of the jet, where it is still very col- 
limated, it can therefore be expected that a large fraction of the protons 
interacts before escaping, producing neutrinos. Farther outside in the jet, 
the size of the knots increases with the distance along the jet and the 
optical depth will decrease. So, even in the case of synchrotron radiation. 
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the main location of neutrino production is expected at the foot of the 
jet, while protons will be able to escape the source at larger distances 
from the core. 

(3) Proton interactions with the proton field 

For FR-1 galaxies with closed tori, the outer edges of the jet will pass 
through the torus, at a distance of Zj 1 — 10 pc from the central black 
hole. The column depth of the torus i^ X ~ 4 ■ 10^^ cm~^. The proton- 
proton optical depth is therefore 

W.. =^-cT„^2-10-3, (11) 



with cTpp ~ 50 mb. This interaction efficiency of 0.2% is small compared 
to the proton-photon optical depths discussed above. In addition, only the 
small fraction of all protons in the outer parts of the jet will interact (see 
Fig. 1). We conclude that proton-proton interactions can be neglected 
herj^. 

To sum up, proton-photon interactions are the dominant source of high-energy 
neutrino production. Optical depths are around unity at the foot of the jet 
and decrease with the distance from the core. Therefore, we expect neutrino 
production in the ffist large shock at around Zj ~ 3000 r^. 



1.3 Consequence for neutrino and Cosmic Ray emission 



As a consequence of the high optical photon-proton optical depth close to the 
foot of FR-1 jets, there has to be abundant neutrino production near that ffist 
shock, and that shock is in the already high speed relativistic flow. On the other 
hand, as we now suspect that there is lateral subs tructure in these relativistic 
flows (iGopal-Krishna et al.l . l2004l: iLovelacd. Il976l. e. g .), th ings could be subtle 
as regards the specific beaming (ILind fc Blandfordl . Il985l . e.g.). This implies 



that the spatial conditions suffice for the production of 10^^ GeV protons at 
the base of the jet, but losses due to proton interactions with the disk's photon 
field will take several powers of ten from that as we discuss in the following sub- 
section. The ultra high energy particles we do obse rve probably receive their 



last energy increase at the last strong shock (see iBiermann fc Strittmatter 



^ For a discussion of the torus' column depth in AGN, see lZier &: BiermannI ( 20021 ) 
and references therein. 

^ This situation could change drastically for AGN with extremely high matter den- 
sities, i.e. GPS/CSS so urces, where both proto n-proton interactions and pion-proton 



O'Dea k Bairn] (Il997l ). Those sources are not considered in 



interactions can occur 
this paper: due to their high proton-proton interaction rate, they are not good can- 
didates for the observed charged Cosmic Ray spectrum. They could be very efficient 
neutrino emitters. 
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( 1l987l )). since the optical depth decreases with the distance from the AGN 
core. Neutrino observations may have the spatial resolution to check on this. 



Disk 



neutrino 
production 

p Y-> ... V 



iog(dN ydE^'E dE^) ~t ~ z 



log(B)~ z 



roton 1 
mission I 




FRS (BL Lac) 



SRS (FR-i) 



-3 -2 -1 +1 +2 +3 



iog(z. /pc) 



Fig. 3. Schematic figure of the AGN jet. The colored regions represent the emission 
regions. The detailed view of the shock structure discussed in Fig. 2 is not shown 
here, but it is assumed implicitly. While neutrino production happens in early, dense 
shocks, protons are more likely to come from the last shock, which is optically thin 
to proton-photon interactions. In this scenario, protons can be observed from FR-I 
galaxies, viewed from the side, and also from BL Lac objects, looking directly into 
the jet. Neutrino emission, on the other hand, is collimated and can only be observed 
from BL Lac objects. 

The schematic view of this AGN model, showing the AGN jet with its shocks, 
is displayed in Fig. 3. Since neutrino production happens close to the foot of 
the AGN jet in a strongly accelerated reference system, the emission of neu- 
trinos is beamed. The neutrino background due to active galactic nuclei has 
been estimated many times, usually in conjunction to the process o f acce l- 
erating ultra high energ y particle s , see e.g. reviews by iGaisser et al.l (119951 ): 
Halzen fc Hoopej ( 2002 ): Becker ( 2008 ). In particular, flat spectrum radio 
quasars (FSRQs), interpreted as FR-II sources with their jet pointing t oward s 



Earth, have been predicted to emit neutrin o s by iMannheim et al.l (119921 ) 



a p. 

Atovan fc Permed mOll l2003h 



Becker et al.l (120051 ). Further, one extensive 



work is by lBednarek fc Protherod (119991 ) . who consider the acceleration of pro- 
tons very close to the disk, up to distances of 10^^ cm, and using a reconnection 
model for the acceleration, a process well studied in activity regions on the Sun. 
In our model, on the other hand, we discuss distances of about 3 ■ 10'^ gravita- 
tional radii, so for typical black hole masses of 10^ Mq, so distances somewhat 
larger, about 4-10^^ cm, and adopt the point of view th at shock wave accelera- 
tion is the dominant process. On the other hand, just as lBednarek fc Protheroe 
( I1999I ). we use the disk radiation as the photon field for interaction. Their fig- 
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ure 3a demonstrates the limit of such an appro a ch: A typical disk tempera ture 
is of order 3 ■ 10^ K (IMalkan &: SargentUl982l: iDonea &: Biermanru. Il996r) . or 
sometimes even higher, as also used in iBednarek fc Protherod (119991 ). They 
confirm that the optical depth is of order unity. However, there is one critical 
difference: In shock acceleration the energetic protons are isotropic in the co- 
moving fram e, here moving with a Lorentz factor of order 10 — 30, or perhaps 



even higher (iGopal-Krishna et al.l . 120071 ) . So we do not have in our concept a 



proton moving rectilinearly straight out, but prot ons in a phase space distr i- 
bution moving along with the relativistic jet flow flFalcke fc Biermannl . Il995l ). 
Hence, particle emission is beamed along the jet axis. 



Neutrinos are therefore only observable from Earth if the jet points (almost) 
directly towards the observer. Since protons are likely to be produced in the 
last strong shock, AGN seen from the side also contribute to the flux of Cosmic 
Rays. 



1.3.1 Flat Spectrum Radio Sources 



As a consequence of the beamed neutrino emission, neutrinos correlated to 
the proton emission of radio galaxies can only be observed from fiat spectrum 
radio sources (FRS), i.e. FR-I galaxies with their jet pointing towards Earth. 

Most recently, the detection of a double-structure flare in optical and X- 
ray wavelengths from the flat spectrum radio source BL Lacertae was in- 
terpreted as the emission from particle acceleration in moving shock fronts 
very close to the central black hole. Then, when the plasma reaches the tur- 
bulent zone which marks the transition from moving to stationary shocks, 
a second sync hrotron flare is s een f rom particle acceleration in the station- 
ary shock, see iMarscher et al.l (120081 ) for details. During the first X-ray flare, 
> 200 GeV emission from BL Lacertae wa s dete cted by the MAGIC tele- 



scope (lAlbert et al.l . 120071 ). iMarscher et al.l (120081 ) explain this high-energy 



component by Inverse Compton scattering of the accelerated electrons with 
the synchrotron photons. An alternative explanation would be the production 
of very high-energy photons via photohadronic interactions. Here, charged and 
neutral pions are produced, the charged pions decaying to produce neutrinos 
and the neutral pions resulting in TeV photons. Thus, if the TeV signal is 
of hadronic nature, i.e. if it arises from tt" —decays, neutri nos are produced 
simultaneously, as also discussed by lMannheim et al.l (Il992l ). 



1.3.2 Steep Spectrum Radio Sources 

Steep spectrum radio sources are AGN seen from the side, often showing 
a detailed jet structure, with radio knots along the jets (FR-I galaxies) or 
radio lobes at the outer end of the jet. We predict that those sources are 
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weak neutrino sources, considering the model presented above: The beamed 
emission from the jet is not directed towards Earth and very strongly focused 
neutrino emission cannot be observed from Earth. 



As for one of the brightest radio galaxies in the sky, Cen A, it will be possible to 
localize the origin of neutrinos within the source: the extension of Cen A is ten 
degrees on the sky, while the resolution of neutrino detection with IceCube will 
be around 1°. So, the flux of high energy neutrinos produced in each subsequent 
shock region will diminish with the distance from the black hole and the 
maximal energy will go up. However, as the jet of Cen A is not pointing near to 
the line of sight, we will see neutrinos only from secondary particles decay, after 
a primary charged particle has scattered in magnetic fields, near the boundary 
or outside the relativistic jet. In this paper, we will calculate the intensity of 
the neutrino signal connected to the possib le correlation of UHECRs and th e 
distribution of AGN as observed by Auger (lAuger Collaboration! . 120071 . 120081 ). 



2 Estimate of the Cosmic Ray flux from nearby sources 



Recent Auger results (lAuger Collaboration! . 120071 . |2008| ) indicate that at least 
20 events above 57 EeV are correla ted with the distribution of n earby AGN in 
the Veron-Cetty & Veron catalog. IVeron-Cetty fc VeronI (120061 ). abbreviated 
VCV catalog in the following. The integrated exposure above 40 EeV is given 
as Exposure = 9 ■ 10^ km^ yr sr and can be assumed to be the same above 
57 EeV. Thus, the integral flux of UHECRs is 



N{> E) 



^{events) 
Exposure 



:i2) 



and in this case 
N{> 57EeV) 



20 



9-103 



km ^yr ^sr ^ = 7 ■ 10 ^^s ^ sr ^ cm ^ 



(13) 



If the flux is correlated to a single or few point source (s) rather than to many 
sources, the field of view (FoV) of Auger l^Augor has to be taken into account 
in combinatio n with the dechnation of the sources uJsnurrpiS), as discussed i n the 



appro aches of lCuoco fc HannestadI (120081 ) ; iHalzen fc O'Murchadhal (120081 ) ; iKoers fc Tinyakov 
(120081 ). As opposed to those models, we calculate the diffuse contribution as 
a conservative estimate, based on a detailed model concerning the physics of 
the AGN jet. 



Based on the correlation claimed by Auger, we assume in the following calcu- 
lation that these 20 events come from AGN in the supergalactic plane (SGP). 
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If the source population for the origin of the calculation is verified to be radio 
galaxies, it is likely that even more than 20 events come from AGN, since the 
VCV catalog is not complete. 



The differential flux from the SGP can be assumed to follow a power-law 
dNcR 

CR SGP 



^SGP ■ -E'en 



(14) 



Here, the cosmic ray energy Equ is given in the laboratory frame at Earth. 
In all following calculations, energies E are in the laboratory frame at Earth, 
energies E^°^'^'^^ = {1 + z) ■ E are given in the laboratory frame at the source 
and E' = i^sourcc^'-p _ £] ■ [1 + z)/r represent the energy in the shock rest 
frame, with F as the bulk Lorenz factor of the shock and z as the cosmological 
redshift of the source. All calculations are done in the way that energies have 
units of [E] = = 1^^'] = GeV. The normalization factor AgQp can 

be determined by using Equ. (13) and comparing it to the integral form of 
Equ. (14): 



N{> EZ 



ger; 



dE, 



^min 

Auger 



dN< 



CR 



CR 



dEcK 

1 / 7-imin 



^SGP / -E'cR dEcK 



Auger 

-p+i 



ISGP-(P-I) (-EXugerJ 

Thus, the normalization factor Agcp can be calculated to be 
^SGP = N{> EZ^) .{p-1). {EZlrr' 



(15) 



(16) 



with ^'Auger in uuits of GeV and [Asgp] = GeV~^s~^ sr"^ cm~^. The spectral 
index p of UHECRs is observed to be 2.7. However, stochastic shock accelera- 
tion itself resu lts in a spectrum with a spec t ral index of 2.3 or e ven flatter, see 
e. g. papers bylBednarz fc Ostrowskil fll998l ): iKardashevI (119621 ): iBaringI (120041 ) 



or 



Meh et al.l (120071 . l2008l ). The steep observed spectrum may be due to the 



fact that there is a distribution of maximum energies or due to propagation 
effects. Therefore, we use index values of p = 2.3 and p = 2.0 for the sources 
of UHECRs in the SGP. For a threshold energy of 57 EeV and an integral flux 
of 7 ■ 10~^^s~^ sr~^ cm~^, the numerical value of the normalization constant is 



A 



SGP 



9 ■ 10^^ GeV-^s^i sr-^ cm'^ for p = 2.3 
4 ■ 10-1° GeV-^s-i sr-i cm'^ for p = 2.0 . 



(17) 



Figure 4 shows the observed cosmic ray spectruni at UHEs, weighted with 
Ep^'"^ . Data are from the HiRes ( Abbasi et al. . 2008bl ) and Auger ( Abraham et al. 



14 



E"'*exp(-E/1e12 GeV) 
E""*exp(-E/1e12 GeV) 
E""*exp(-E/1el2 GeV) 




8 8.5 9 9.5 10 10.5 11 11.5 12 

loglO(E/GeV) 



Fig. 4. Differenti al UHECR spec t rum w eighted with E 



-2.7 



measurements from 



Auger (circles ), ( Abraham et al. . 20081 ). and HiRes-I and HiRes-II (triangles), 



(jAbbasi et al.l . l2008bl ) . The data from the two different experimen ts are comp arable 



when considering small systematics in the energy measur;ments «H). The 
differential flux from the supergalactic plane calculated for E~'^ (black line), E~'^'^ 
(red line) and E~'^''^ (blue line). An exponential cutoff was applied at 10^^ QgV. 



20081 ) experiments. The flux from the supergalactic plane as calculated above 
is indicated as the blue {p = 2.7), red {p = 2.3) and black {p = 2.0) lines. The 
three indices represent different models for the description of Cosmic Ray pro- 
ductio n and propagation following the theor etical results of iBerezinskv et al. 



( 2006 ) in the case of p = 2.7, the work of ( Bedriarz fc Ostrowskil. 



(120071 . 



1998 ) for 



20081 ) for 



p = 2.3, using parallel shocks, and the results of iMeli et al. 
oblique shock conditions [p = 2). The energy flux can be calculated from the 
differential spectrum as 



thJ 



EdN/dEdE. 



The total energy flux measured in cosmic rays above -Eth = 10^'^ GeV is 
approximately 

jtot-CR ~ 10-' GeVs-^ sr-i cm'^ . (19) 



The corresponding energy flux from the supergalactic plane can be calculated 
by using the same energy threshold and the differential spectrum determined 
above: 
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JSGP~CR ■ 



^ ■ for p ^ 2 

^SGP ■ ln(-E'max/-E'th) 

4 ■ 10-9 sr-i cm-^ for p = 2.3 

9 ■ 10-1° GeVs-i sr-i cm-^ for p = 2.0 



(20) 
(21) 



Here, a maximum energy of -Emax = lO^^'^ GeV was used, taking into account 
the GZK cutoff. This imphes that about 1% of the total cosmic ray flux above 
the ankle is made up by sources in the supergalactic plane. 



3 Neutrino flux estimate from Auger measurements 



With protons being accelerated in AGN iets, ne utrinos can be produced in 
photohadronic interactions (see e.g. iBeckerl (120081 ) for a review). 



PI 



A' 



TT^ p 



n 



7 7(p) 



2/3 of the cases 



/i+ (n) Ue '^fi (e~'' n) 1/3 of the cases 



.(22) 



Here, the branching ratio for charged pion production is 1/3. About 1/2 of the 
pion's energy goes into the sum of muon and anti-muon neutrinos (Notation: 
u := + V^). In addition, neutrinos oscillate on their way to Earth from a 
ratio of 

(z/e, v^,, i^r)|sourco = (1 : 2 : 0) to (i/^, v^, z^r)|Earth = (1:1:1), (23) 



see e.g. IStanevI (120041 ) and references therein. Given an optical depth Tp^ for the 



production of the Delta- resonance, the total neutrino energy flux jy, radiated 
in the solid angle VL^, {VLy Jl) = I El dNyj dE^ dE^, and the total proton energy 
flux, focused within the sohd angle flcR, (^CRjp) = / EpdNp/dEpdEp, both 
given in the shock rest frame, are therefore connected as 



{Qy ■j'u)=Tp^-^-yy i^CR fp) 



P7 
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(24) 



The optical depth will be set to = 1 in the following calculations (see 
Section 1 for a detailed discussion). This correlation can be used to estimate 
the neutrino flux to be expected from the supergalactic plane, given the proton 
flux at Earth as in Equations (14) and (16). First, we need to connect the 
energy flux in the shock rest frame and the one measured at Earth. 
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3. 1 The energy flux 



The diffuse energy ffux of a certain particle species is given in terms of tlie 
energy spectrum at Eartli as described in Equ. (18). In terms of the single 
source spectrum, it can be written as 



''"'^^ dN 1 d'^n dV 



^max -LTnax 



3= / / dE'"''""'dzdLE'°''''''— ^-JT- 377-77 ■-^•(25) 



_g;sourco Zn-,in L„ 



Here, we are taking into account the following facts: 

(a) Each source contributes with dA^/dE^""''^'' 

(b) The flux of each source decreases with Aird'j^, where di is the source's lu- 
minosity distance. 

(c) The source number per comoving volume dV/dz and per luminosity inter- 
val, dn/dV/dL - for protons, we use the radio luminosity function of FR-I 
galaxies, observed neutrinos are only produced in FRS. 

(d) Sources up to a redshift ^max, with an absolute upper limit of 0.03 as the 
outskirts of the supergalactic plane contribute to the energy flux, both for 
neutrino and for proton sources. The minimum redshift is given by the 
distance of the closest plausible source, Cen A, at = 0.0008 in the case 
of FR-I galaxie£3 (CR sources), and Perseus A at z^^^ = 0.018 in the case 
of FRS (Neutrino sources). 

(e) The luminosity integration limits for AGN have been chosen to match the 
observed distribution of FR-I galaxies, L^im = 10^° erg/s and Lmax = 
]^q44 QYg/s. This applies for FRS as well, as these are generally believed 
to be a sub-class of FR-I galaxies. 

As we would like to compare the diffuse energy flux at Earth to the energy 
flux in the shock rest frame, we use iJ'^ource = y ■ E' and receive 



r 

Air 

r 

An 



dE' E' 



■3 -n, 



dN 
dE' 



•Smax irnax 



dz dL 



d'^n dV 



d\ dV dL dz 



(26) 
(27) 



Cen A is located at a distance of 3.5 Mpc which corresponds to a corrected 



redshift of 0.00083 for a cosmology of VLrn = 0.3, f^A = 0- 7 and h[ 



approximate values from most recent WMAP-5year results (jKomatsu et al 



0.7 , usin g 
20081 1. 



The apparent, measured redshift is much higher than the actual cosmological one, 
due to the gravitational influence of the Great Attractor on nearby sources. 
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with 

J J Airdl dVdL dz 

2min Lmin 



Inserting Equ. (27) into Equ. (24) for both neutrinos and protons, we have 

1 ficR . 

J^ = T77T^ F; -^CR- (29 

12 TcR 12,, ncR 



3.2 Discussion of parameters 

• Cosmic ray energy flux at Earth 

The proton energy flux can be written as 

isGP-CR = ^SGP / dEcKEX (30) 

f(p-2)-i-(i?cR""°)"^+' forj9^2 
= ^'"^^"|ln(^) forp = 2 ^''^ 

= ^ ■ ^Scr) ■ (^Ser)""' " (^CR )""^' ioT p 2 .(32) 

The proton energy is given in units of GeV and [Ap] = GeV~^s~^ sr~^ cm~^. 
For protons, the lower energy threshold is given as the proton mass boosted 
by TcR, E^^ = ^ TcR ■ 1 GeV. 



Neutrino energy flux at Earth 

It is assumed that the neutrino spectrum traces the proton spectrum, 
dN^ 



dE„ 



A. ■ E;^^ (33) 



with ai, ~ p, El, in units of GeV and [A^] = GeV ^s ^ sr ^ cm ^. The 
neutrino energy flux can then be expressed as 

j, = A, J dE,E;''^+' (34) 

11 

- 2)-l ■ (E^m)-a,+2 ^ . r-a.+2 fo^ 2 



(35) 

In ( ) for a,, = 2 . 



A„ < 



The lower energy threshold of one fourth of the pion mass, boosted by 
Tiy, E'™™ = ■ m^/4 = Fy ■ 0.035 GeV, dominates and the upper energy 
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threshold can be neglected in the case oi ^ 2. For = 2^ the logarithm 
for neutrino energy equals approximately the one for proton energies and 
the factors cancel in the end. 

Concerning the threshold energies a nd spectral behavior of the flux, our 



model differs from other approach es by lCuoco fc HannestadI (120081 ) ; iHalzen fc O'Murchadha 
( 120081 ) ; iKoers fc Tinyakovl (120081 ) in the sense that we assume proton inter- 
actions with the photon field f rom the disk and from the synchrotron radia- 

tion in t he jet. In the model s oflCuoco fc HannestadI (120081 ) ; iHalzen fc O'Murchadha 



( 120081 ): iKoers fc Tinyakovl (120081 ). X-ray photons are assumed to interact 



with the protons in the sources, which results in a broken power-law behav- 
ior when calculating the neutrino flux. The reason is that the optical depth 
for proton-photon interaction changes with energy, and this change happens 
at a neutrino energy of around ~ 10^ GeV for X-ray photons. In the case of 
optical or radio photons, the break is at very low energies, so that it is not 
relevant for our calculations. Thus, the neutrino flux calculated here follows 
a single power-law just as the proton spectrum does. 



Lorentz factor 

As neutrinos come from early shocks near the black hole, the beaming factor 
is typically stronger than for protons from late, outer shocks. We will assume 
in the following that the boost factor for neutrino production is a factor of 
~ 3 higher than for protons, 

"T 

3. (36) 



CR 



Solid angle 

We assume in the following that Q cc 6"^ (x 1 /F^ as the typical opening angle 
for relativistic sources. In particular, Q„ ^ and f2cR oc 1/F^j^. 



Redshift dependence 

The redshift dependence of the neutrino energy flux is given as 



z" T 

max ^ 



max -^^max 



Uy = I I dzdL ^ ^ 



indl dVdL 

20.018 Lmin 



FRS 



dz 



(37) 



In a first approach, we will calculate the total diffuse neutrino flux from the 
supe rgalactic p l ane a nd use z^^-^ = 0.03 as the upper integration limit, see 



pl ane a n( 

also iDas et al ] (l2008h . The radio luminosit y function for fiat spectru m ra- 
dio sources, d'^n/dV/dL\p-^^, determined by Dunlop &: Peacock ( 1990l ). The 



authors do not distinguish between the high- and low-luminosity sources. 
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so both FR-II type objects ("FSRQs") and FR-I ("BL Lacs") are included. 
Since the low-luminosity part of the sources makes up most of the popula- 
tion, we neglect the contribution of FSRQs here and interpret this function 
of fiat spectrum radio sources (FRS) as an approximation of the FR-I com- 
ponent of blazars. 

As discussed before, neutrino emission is beamed and originates from FRS 
only. For protons, all FR-I galaxies contribute: 



dz dL 



1 



d'^n 



0.0008 L„ 



471 dl 



dV dL 



FR-I 



dV 

dz 



(3J 



The r adio luminosity function of FR-I galaxies is given in IWillott et al 
(120011 ). As the Auger correlation at most reveals a general correlation be- 
tween the distribution of AGN and cosmic rays, the exact number of sources 
is not known. Note that the extension of sources goes with the distance 
squared rath er than with the total volu me due to the disk-like structure of 



the SGP, see ICavaliere fc Mencil (119971 ). This is compensated by using the 



volume integration for both the neutrino and the Cosmic Ray sources: The 
relevant ratio n^/ncK gives the correct result. 

Figure 5 shows how the ratio rij^/ncR depends on the maximum redshift 
of contributing sources to the observed correlation. The more sources con- 
tribute, the larger ncR and the smaller the ratio. The most conservative as- 
sumption is to assume that all sources in the supergalactic plane contribute 



to the correlation, = 0.03, see also iDas et all (120081 ). The most opti- 



mistic assumption is that only the nearest sources contribute, z^^ = 0.001. 
Thus, the values lie between 



0.1 < 



< 5. 



'^CR 



(39) 



3.3 The diffuse neutrino flux from the supergalactic plane 



Inserting the results from Section 3.2 into Equ. (29) yields a numerical value 
for the neutrino normalization constant, 



A, 



5— p 



'^CR VrcR/ 12 

777 

N{> E^) .{p-1). i-^y-' ■ [E^^J^-^ , (40) 
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Fig. 5. Redshift factor n^/ncR for with variable maximum redshift for ncR- If only 
very nearby sources are responsible for the correlation of UHECRs with the dis- 
tribution of AGN, small as zz^^ = 0.001. If the correlation comes 
from the entire distribution of sources, the maximum redshift is given by the exten- 
sion of the supergalactic plane, ^^^ax — 0.03. The minimum redshift for integration 
is taken to be z'^^ = 0.0008 for the cosmic ray factor, using Cen A as the clos- 
est contributing source, and the neutrino factor has z^;^ = 0.018 as Perseus A as 
the closest contributing source. The maximum redshift for neutrinos is taken to be 
^max — 0-03. 

where r^^ = 1 was assumed, as discussed in more detail in Section 1. 

The neutrino flux for an E~'^— and an E'"^'^— shaped spectrum is shown 
in Fig. 6, assuming that the Auger correlation is caused by all sources up 
to z^ax ~ 0.03. The neutrino flux is higher at lower energies for the steeper 
spectrum E^'"^'^. The reason is that the normalization is done at Auger's 
threshold energy. If the spectrum is steep ~ E^~'^'^^~'^'''' , the energy content at 
lower energies is higher. If the spectrum is very flat (e.g. E~^), the contribution 
at energies lower than Eq^ is also lower (see Fig. 7). If the spectrum flattens 
towards smaller energies, this reduces the normalization as well. 

The numerical value of the neutrino normalization constant for an E'"^— spectrum 



IS 




0.03 



0.002. 



(41) 
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The expected neutrino flux is about a factor of 35 higher if the closest AGN 
produce the correlation of cosmic rays and AGN. 




Fig. 6. Diffuse neutrino flux from the supergalactic plane, based on Auger data. 
Tlie black line, labeled ®, is based on an spectrum , the red line, labele d ®, 
is calculated using E~'^'^ . AMANDA data are taken from Miinich et al. I (|2007l l and 



the AMANDA limit is given in lAchterberg et al.l \2mi \. 
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Fig. 7. Proton energy density, based on an ^ spectrum (black line). Assuming 
that a spectrum steeper than at the highest energies actually approaches 

at low energies, the energy density saturates at p = 2.3 (red, dashed line). 
Still, there is a difference in the energy density of two orders of magnitude between 
assuming an E~'^ spectrum and an E~'^'^ spectrum. 
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3-4 The total, extragalactic diffuse neutrino flux 



The total, extragalactic neutrino flux can be calculated by integrating up the 
neutrino flux redshift dependence u p to z';!^^^ = 7, where the first a ct ive galax- 



i es ar e beheved to contribute, see iBouwens fc lUingworthl (120061 ): live et al. 



( 120061 ) ■ The result is shown in Fig. 8. Here, we use the conservative assump- 
tion that sources up to = 0.03 contribute to the correlation between cos- 
mic rays and the distribution of AGN. For an shaped spectrum (black, 
horizontal, solid line), the fiux is a bout a factor of 10 higher than the pre- 
diction made in iBecker et al.l (120051 ) . where the jet-disk symbiosis model was 
used to estimate the neutrino fiux from FSRQs. Since we normalize the AGN 
spectrum at the highest energies, a steeper neutrino spectrum (e.g. E~'^-^, red, 
solid line) leads to a higher fiux at low energies. Although the figure seems to 
indicate that an E~^'^— fiux is already excluded by AMANDA data, it should 
be noted that the limit is derived for an E~'^— shaped spectrum. To be able 
to exclude the spectrum, the limit needs to be calculated fo r the same 



spectral behavior, see (IBecker et al.l . l2006l : lAchterberg et al.l . l2007l ) 




Fig. 8. Diffuse neutrino flux from flat spectrum radio sources up to z < 7. The black 
line, labeled ®, is based on an £'~^ spectrum, ttie red line, labeled @, is ca lculated 
using ^-2. 3 'YYiQ lines label e d @ and ® are shown for comparison, from ISteckei 
and iMannheim et alJ I (|2nnih . respectively. The two models use the high-en- 
ergy component from AGN, > MeV, to ca, l culate the correlated neutrino flux. 
AMAN DA data are taken from iMiinich et aD (j2007l ) and the AMANDA limit is 



Achterberg et al.l (|2007l '). The dashed lin e, labeled ®, shows the prediction 



given m 

of the neutrino flux from FSRQs as calculated in IBecker et al.l (j2005l ). 
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3.5 Uncertainties in the determination of the spectrum 



Neutrino flux calculations typically bear three significant sources of uncertain- 
ties, the first one being in the total normalization of the spectrum, the second 
one lying in the uncertainty of the spectral behavior. The third one is the 
maximum neutrino energy of the source class. Both quantities always rely on 
the internal properties of the source class, which are typically poorly deter- 
mined. This is true in general and does not only apply to these calculations. 
Here, we explain how these calculations are effected by the uncertainties. 



(1) The normalization 

The uncertainty in the normalization of neutrino spectra has three main 
components: the measured spectrum used to normalize the neutrino flux, 
the optical depth in the sources and the opening angle of neutrino and 
proton emission. 

The correlation of the highest energy Cosmic Rays with the VCV catalog 
is based on 20 events. These statistics need to be enhanced in order to 
achieve a more precise prediction of the actual Cosmic Ray flux from the 
source class. In addition, the composition of the spectrum is important 
for the neutrino flux, since protons produce more neutrinos than heavy 
nuclei do. 

Secondly, the optical depth of the sources depends on the size of the accel- 
eration region, on the luminosity, on the boost factor and on the photon 
density in the acceleration region. Those properties are known for a few 
single objects, but they can vary with the source and also with time in a 
given source. 

The ratio of opening angles of proton and neutrino emission is conser- 
vatively taken to be ~ 3. However, the real ratio clearly depends on the 
properties of the single AGN, on the location of the first and last shocks. 

(2) The spectral index 

Depending on the orientation of the shock towards the magnetic field and 
the boost factor, the spectral behavior can vary. For parallel shocks, a 
spectr al behavior of E'"^'^ is expected as discussed by lBednarz fc Ostrowski 
(119981 ). If using large angle scattering instead of pitch angle scatter- 
ing for parall e l sho cks, particle spectra of up to E~^'^ can be obtained 



Stecker et all (120071 ). For oblique, subluminal shocks, the spec tra behave 



as E ^-^^ depending on the boost factor, as discussed by lMeli et al. 
(l2008h . 
(3) Maximum energy 

As discussed before, the maximum energy of the protons, and hence of the 
neutrinos, depends on intrinsic quantities like the magnetic field, boost 
factor and the disk luminosity. 
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So, to conclude, the importance of large volume neutrino detectors is enhanced 
by the fact that protons or high-energy photons cannot give unambiguous 
evidence for the spectral index, nor for the strength of the neutrino flux. In 
order to identify the nature of the shocks in not only AGN, but any Galactic or 
extragalactic accelerator, neutrinos are essential. In this paper, we predict the 
region of neutrino production in AGN jets. We predict that fiat spectrum radio 
sources should be dominant. This can easily be tested by future experiments 
like IceCube and KmSNeT. 



4 Summary and implications 



In this paper, we present a model for Cosmic Ray and neutrino emission 
from active galactic nuclei. A first evidence for the correlation of the observed 
UHECRs above 60 EeV and the distribution of nearby AGN gives rise to 
the prediction that UHECRs may come from FR-I galaxies. Although FR-II 
galaxies would be a good candidate due to their high intrinsic luminosity, they 
are not abundant enough and too far away to lead to such a correlation. FR-I 
galaxies, on the other hand, have a much higher source density and there exist 
several nearby sources which could be responsible for a directional correlation. 
The most prominent FR-I galaxies are Cen A and M 87, as well as the fiat 
spectrum radio sources BL Lac and Perseus A. Flat spectrum radio sources 
with low luminosities, called BL Lac objects, are typically interpreted as FR-I 
galaxies with their jet pointing directly towards Earth. 



Observations of the je t structure of M 87 (I Walker et al.l . l2008l ) and BL Lac 
( iMarscher et al.l . |2008| ) near the central black hole indicate that, while there 
can be moving shocks between 10 and 1000 Schwarzschild-radii, the first 
strong, stationary shock occurs a t ^ 3000 Schwarzschild-radii 



^ , as al- 



ready discuss e d by Markoff et al.l (120011 . |2005| ) . This is further confirmed by 
Britzen et all (|2008|) in the case of the BL Lac type object S5 1803+784. At 



~ 3000 Vg, the optical depth for proton interactions with photons from the disk 
is around 2%. With synchrotron photons from the jet as a target, the optical 
depth is ~ 90% at the same distance from the central black hole. In both cases, 
the optical depth decreases with the distance from the black hole, Tp^ ~ z~^. 
Therefore, we predict that neutrinos are produced in the narrow jet, close to 
the central black hole. Their emission is beamed due to the proton's highly 
relativistic motion along the jet. Therefore, neutrinos can only be observed 
for sources pointing directly towards Earth. Protons, on the other hand, dom- 
inantly arise in the last strong shock of the jet, at several kpc distance from 
the central black hole, where the optical depth is close to zero. 



Using the correlation between UHECRs and the distribution of AGN as a 
measure for the Cosmic Ray flux from the supergalactic plane and connecting 
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this to the neutrino emission produced near the foot the AGN jet results in 
a predicted neut rino signal about an ord er of magnitude below the current 
AMANDA limit lAchterberg et all (j2007l ). IceCube will be able to test this 
model within the first years of operation. 

The general terms of this model are independent of Auger data, and cou l d hav e 
been normalized to the established UHECR flux (IGaisser fc Stanevl. 120061) . 
assuming that radio galaxies ar e the sources (iGinzburg fc Syrovatskiil . Il964l : 
Biermann fc Strittmatterl . 119871 . e.g.). Au ger does conflrm a correlation wit h 
the distribut ion of active galactic nuclei (lAuger Collaboration! . 120071 . l2008l ). 
while HiRes ( lAbbasi et all . l2008al ) does not, using the same energy threshold, 
and the same procedure. It is clear that more data are required, and the 
energy threshold might ultimately be seriously different; the results show that 
a conflrmation from a larger data set will be necessary. 
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